Simulating an interacting gauge theory with ultracold Bose gases 
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We show how density dependent gauge potentials can be induced in dilute gases of ultracold atoms 
using light-matter interactions. We study the effect of the resulting interacting gauge theory and 
show how it gives rise to novel topological states in the ultracold gas. We find in particular that the 
onset of persistent currents in a ring geometry is governed by a critical number of particles. The 
density-dependent gauge potential is also found to support chiral solitons in a quasi-one-dimensional 
ultracold Bose gas. 
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Introduction. Our understanding of the fundamental 
interactions between elementary particles is founded on 
gauge fields. The role of the gauge field is to mediate 
the interaction between particles. The simplest exam- 
ple we know of is electromagnetism where charged par- 
ticles interact through exchanging virtual photons. The 
Coulomb potential between charged particles is encom- 
passed by this gauge theory which can be recast in the 
familiar form of Maxwell's equations in the classical limit. 
For this to hold, the gauge field must be dynamical. 
In other words we must be allowed to construct a La- 
grangian which describes the propagation of the gauge 
field in vacuum. 

Gauge theories are not restricted to electromagnetism 
only. The interactions in nuclei are governed by more 
complicated objects as far as gauge fields are concerned. 
There one has to use higher dimensions which typically 
requires a non-Abelian theory, such as the Yang-Mills 
field for the gluons. The complexity of the system in- 
creases dramatically. One of the main goals here is to 
study quantum chromodynamics, which is believed to 
explain the confinement of quarks and chiral symmetry 
breaking at high temperatures and densities. Experimen- 
tally the phase diagram of QCD cannot be mapped out 
completely due to the energy scales and the densities re- 
quired. There are no accelerators powerful enough for 
that. More worrying perhaps, the phase diagram is also 
intractable from a theoretical and calculational point of 
view due to the complexity and the fact that QCD, which 
is highly non-perturbative p] , belongs to the class of ex- 
ponentially hard problems [2]. The solution could be to 
design a special purpose quantum simulator [3] . For this 
we ideally need to be able to create gauge fields which 
are dynamical. 

Very recently, the first few theoretical proposals in this 
direction have appeared @H7], where it was shown that 
it is in principle possible to simulate a dynamical gauge 
theory using cold atoms trapped in optical lattices. Also 
smaller steps towards the ambitious goal of simulating 



QCD using possibly less demanding experimental tech- 
niques may provide some important insights (see for in- 
stance [SHlOj). A more modest problem that generated 
intense interest in the late 1990s was the quest for find- 
ing a pure gauge theory with solutions given by the one- 
dimensional analog of the well-known two-dimensional 
anyons [11] . The first attempt in this direction 12J failed 
to describe one-dimensional anyon solutions [13] , but the 
associated semiclassical, non-linear model of the inter- 
acting gauge theory supported chiral solitons, as shown 
by Aglietti, Griguolo, Jackiw, Pi and Scminara (AGJPS) 
in Ref. [14]. The generation of chiral solitons is clearly 
also an interesting goal to pursue in its own right due 
to the unconventional coherent transport mechanisms in 
the superfluid regime. 

In this Letter, we show how properly engineered laser 
fields can induce an effective density-dependent vector 
potential in a weakly-interacting ultracold Bose gas, 
which constitutes the semiclassical limit of an interacting 
gauge theory for bosons. When the system is tightly con- 
fined such that it forms a quasi-one-dimensional gas, it is 
described, in a one-to-one fashion, by the AGJPS gauge 
theory 14J . We illustrate the effect of such a gauge field 
by studying the one-dimensional dynamics of the free ex- 
pansion, the onset of persistent currents in a ring, and 
chiral solitons in a Bosc-Einstcin condensate. 

An effective interacting gauge theory. There are a num- 
ber of ways to induce artificial magnetic fields in ultra- 
cold atomic gases. The most straight forward method 
which has been used extensively for vortex creation in 
Bose-Einstein condensates relies on stirring the cloud by 
a laser spoon or using asymmetric external traps |15) . 
In the rotating frame the cloud experiences a homoge- 
neous artificial external magnetic field which is propor- 
tional to the rotation frequency. Alternatively one can 
use laser assisted tunneling in optical lattices which in- 
duce the required phases for the tunneling amplitudes 
between the different lattice sites [TBI Hi]- F° r ultra- 
cold atoms optically induced gauge potentials can also 
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be created based on dark state dynamics [T5H20] or Ra- 
man transitions [2T1 [22] . These gauge potentials all have 
in common that they are static and given by the exter- 
nal rotation frequency or laser parameters; there is no 
dependence on the density of the atomic cloud in the 
gauge potential using these techniques. Here, we show 
how a density-dependent vector potential can arise in 
a weakly-interacting Bose-Einstein condensate based on 
collisionally induced meanfield shifts of the electronic lev- 
els in the atoms, which also constitutes the semiclassical 
limit of an interacting bosonic gauge theory. This can 
be done by considering a gas of optically addressed two- 
level atoms forming a Bose-Einstein condensate (BEC) 
with internal state space given by |1), |2). Alkali atoms 
are usually good candidates for this, but fast spontaneous 
decay from these states might render them unusable. On 
the other hand, a setup relying on dark states [33J could 
circumvent this problem. Hence, good candidates for ex- 
perimentally realizing the effects discussed here would be 
for instance strontium, calcium or ytterbium, which can 
have radiative life times of several seconds |24) . 
The microscopic TV-body Hamiltonian is given by 
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is the Hamiltonian for the light-matter interaction, Vi is a 
single-particle external potential, Ik,... is the identity op- 
erator acting on the subspace excluding particles k, . . ., 
and V M = di&g{g n ,gi2,gi2,92-z]S(rk-ri) is the two-body 
interaction with coupling strengths gty = 4Trh 2 a,ij /m. In 
Eq. fl is the two-photon Rabi frequency character- 
izing the light-matter coupling, fa = 0(r/-) is the laser 
phase at particle fc's position, and the laser detuning from 
the atomic resonance is chosen to be zero for simplic- 
ity. However, the meanfield terms stemming from Vk,i 
will introduce an effective detuning. The corresponding 
Hamiltonian which takes into account collisional mean- 
field effects is then given by Hop = p 2 /2m ® 1 + V + U, 
where V = (l/2)diag[ffupi + 312^2,522^2 + guPi], with 
Pi = \^i\ 2 (i — 1)2) the density of population in the 
atomic state i, such that (H)q, GP — {^\Hgp\^)^ where 

l*GP>=®£i|**>. 

For weakly interacting atoms, the coupling strength 
Ml is typically much larger than the mean-field energies. 
The zero-th order approximation to the state of the sys- 
tem is chosen as the usual starting point in Bogoliubov's 
theory of the Bose gas. In this regime, to diagonalize 
U + V we treat V as a small perturbation to U. If we 
define the densities in the dressed states as p± = l^il 2 , 
the corresponding eigenstates of U + V are given by the 
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with eigenvalues gp± ± Ml/2, g = (311 + g 2 2 + 2gi 2 )/4 
and |x± 0) ) = (|1) ± e^\2))/V2, together with the adi- 
abatic approximation such that either p~(v, t) w or 
p+(r, t) sa 0. A general state can consequently be writ- 
ten like \<p) = Ei={+,-} * 2 ( r ^)IXi)- By projecting onto 
one of the dressed states, \x±), we obtain the effective 
Hamiltonian [231 
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where W = ^|(x-|Vx+)| 2 is a scalar potential and 
A± = ih{x±\Vx±) is a geometric vector potential that 
arises from the projection of the full system onto one of 
the dressed states. In order for the adiabatic approxima- 
tion to hold we must ensure that any induced detuning is 
small compared to the Rabi frequency J7. The resulting 
vector potential is then given, to leading order, by 



A± = A(°)±a lP± (r) 



(5) 



where A^ ) = — fV<^> is the single particle contribution 
and ai = (V(j>)(gii — <?22)/(8^) controls the effective 
strength of the density-dependent vector potential. 

In order to derive a meanfield Gross-Pitaevskii type 
equation we apply the variational principle 5C/5^* = 
to the action C = (*f?\(iHdt — H±)\^f), with respect to 
This is in contrast to the single particle situation 
[23] where a Schrodinger equation is obtained after the 
projection onto one of the dressed states. We consider in 
the following the + branch in H± without loss of gener- 
ality and consequently drop the ± index in p±, ty± and 
A± . The resulting equation of motion is then 



(P-A) 

2m 



ai • j + W + gp 



(6) 



where A is given by Eq. ^ together with a nonlinearity 
in the form of a current, 
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The meanfield scalar potential W is given to leading or- 
der by W = |A(°)| 2 /2m. 

One- dimensional physics. The density-dependent vec- 
tor potential gives rise to a number of interesting and 
counterintuitive scenarios. To illustrate this we will in 
the following assume that the cloud of atoms is tightly 
confined such that any motion in the transversal direc- 
tion is frozen out and the dynamics is well described by 
an effectively one dimensional meanfield description. We 
choose (j> = kx as the phase of the incident laser, together 



3 



with the transformation ty(x) — e lkx / 2 %j;(x), which re- 
sults in the equation 
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where W = h 2 k 2 /8m, and a\ — k(gu — g22)/8ClSt char- 
acterises the strength of the current nonlincarity. The 
effective transversal area of the ID cloud is given by S t - 

The aij term can be made influential provided that 
the meanficld shift is relatively large. In the following 
we will illustrate three scenarios, which can be achieved 
by for instance the combination of an average density 
of 4.0 x 10 14 cm -3 , the difference in scattering lengths 
an — 0-22 = 5nm using for instance optical Feshbach res- 
onances [25H28] . and the ratio Ml /En = 10 where /-.';,■ is 
the recoil energy. This results in the dimensionless pa- 
rameter Nai/h = 12 where N is the number of particles, 
which captures the strength of the current nonlinearity. 

Ring geometry. We envisage a tightly confined cloud 
in a ring trap positioned in the x-y plane such that any 
transversal motion can be neglected. In addition we use 
a laser beam propagating in the z-direction which carries 
an orbital angular momentum with ip = £8 where I is an 
integer. This configuration gives rise to a gauge potential 
in the azimuthal ^-direction, hence the situation is simi- 
lar to the linear ID case, but now with periodic boundary 
conditions. The time-independent Gross-Pitaevskii equa- 
tion on the ring of radius R is obtained from equation 
^ by setting x = R9 and ip(x,t) = ^(x) exp[— iEt/H]. 

The solutions are given by ip(9) — \J ^j^e lq6 with nor- 
malisation condition R f ™ d9\£(0)\ 2 = N, where N is 
the number of particles in the ring. The energy differ- 
ence between two different angular momentum states can 
readily be calculated, 
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where q is an integer number which labels the quan- 
tized rotation of the cloud. We see from Eq. ([9| that 
the ground state configuration becomes a function of the 
number of particles. At a certain critical density 
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the ground state changes from one rotational state to an- 
other with q — > q+1. This is in contrast to the standard 
situation where a ring condensate is subject to a static 
homogeneous magnetic field with the onset of rotation 
given by the rotation frequency or laser parameters. 

The free expansion. A numerical solution of Eq. Q 
shows that the free expansion is no longer symmetric 
(see Fig. [lj. In addition the current term induces a 
drift which is proportional to a\ times the density of 
the BEC. The onset of a drift can be explained using 
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FIG. 1: (Color online) Snapshots of the free expansion of a 
harmonically trapped BEC with (gN/ S t )(2m/h 2 k) = 100 and 
trap frequency uj t = hk 2 /m. The expansion is asymmetric 
where a change of sign in 01 changes the direction of the drift. 
The coupling strength for the gauge field was Nai/h = 12. 
Length is in units of 1/fc and time in units of 2m /hk 2 . 



as variational ansatz the solution of a freely expanding 
wavepacket where we allow for a drift velocity £$ of the 
centre of mass, 
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The spatially varying phase is given by S = mxo(x — 
xo)/h, a x (t) = o-Qy/l+ t 2 /t 2 with r = 2m/k 2 h is the 
time dependent width of the gaussian and N is the num- 
ber of particles. From Eq. (11) and the corresponding 



Lagrangian density we obtain an equation of motion for 
the position Xo(t) of the wave packet, 
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Equation ( 12 ) and its solution 
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gives us a way to understand the effect of the dynamical 
gauge potential on the condensate as it expands. The 
increasing width as a function of time drives the drift of 
the center of mass coordinate rro- 

Figure [l] shows that the presence of the current term in 
^ causes the free expansion of the condensate to experi- 
ence a drift that depends on both the sign and magnitude 
of the strength of the dynamical gauge field captured by 
the parameter a%. The onset of a drift can also be under- 
stood as an effect of the asymmetric coupling of the dif- 
ferent momentum components in the initial wave packet 
with the density of the cloud. 

Chiral solitons. Our semiclassical gauge theory, whose 
equation of motion is given by Eq. pi), supports chiral 
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solitons. Our model is equivalent to the AGJPS model 
[14j . with the additional non-linear interaction term gp. 
The existence of chiral soliton solutions is also ultimately 
a consequence of the breakdown of Galilean relativity 
in the corresponding microscopic version of our gauge 
theory. 

We begin by using the gauge transformation, 



2.5 
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dx'p{a/,t)-iWt/h)$(x,t). (14) 



Eq. ([8]) then simplifies to 
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where the gauge-transformed current becomes 
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Equation ( 15 ) can be solved by first writing the wave 



function in the form [29] 

$(x, t) = £(x - ut y(urnx-(\ m u 2 +nt))/h^ ^ 

where £(x — ut) is a real valued function and p is the 
chemical potential. The current consequently transforms 
into j(x) = ut, 2 where u is the speed of the soliton. The 
resulting differential equation for £(x — ut) is 



(18) 



For g — g — 2a\u > we find in particular the dark 
soliton solution 



$(x, t) exp[i(mu(x — ut) — (l/2mu 2 + p)t)/h] 
coth[(x - ut)/(y/2l )] 



(19) 



where po is the background density, p = gpo and Iq = 
hj ^J2mgpQ. For g < we obtain a bright soliton 



$(x,i) exp[i(mu(x — ut) — (^mu 2 — p)t)/h] 
cosh((x - ut)/(y/2l )} 



(20) 



with p = \g\po/2. The solutions in equation (191 and 



(20) are chiral, which means that the solitons can only 
propagate in a specific direction for a chosen velocity. 
Interestingly, if g = 2a\u we are in a situation where 
the current non-linearity cancels the mean field interac- 
tions between particles, with no soliton solutions present. 
Depending on the precise physical setup this particular 
situation may or may not be possible to reach due to 
a breakdown of the adiabatic assumption or a violation 
of the perturbative assumption. The concept of a chiral 
soliton can be illustrated by considering the reflection of 
a BEC from a hard wall. In Fig. [2] we show how a bright 
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FIG. 2: (Color online) The normalised width a{t)/a{0) of 
the bright soliton. The blue solid line shows a(t)/a(Q) with 
current strength Nai/H = 4.0 and (gN/ St)(2m/h 2 k) = 0.5, 
whilst the red dashed line indicates the width of an initially 
identical soliton without the current nonlinearity (Nai/h = 
0). After reflection the soliton starts expanding due to the 
change in nonlinear strength. The inset shows snapshots of 
the density of the soliton prior to (solid blue) and after reflec- 
tion (dashed black) at times t — 0,2.0,3.5,5.0. All lengths 
are in units of 1/k and time in units of 2m/ hk 2 . 



soliton initially moving in the positive x-direction is de- 
stroyed after reflection. A standard bright soliton would 
retain its width a(t) = *J (x 2 ) after reflection whereas the 
chiral soliton is found to start to expand after reflection. 
The change in the nonlinear strength due to the change 
in momentum after the reflection results in a state which 
is not the soliton solution any more, hence the solution 
is no longer confined. 

Conclusions. In this Letter we have shown how an in- 
teracting gauge theory for a BEC can be generated. The 
equation of motion includes a current non-linearity and 
in the quasi-one-dimcnsional regime our model is iden- 
tical to the Aglietti-Griguolo-Jackiw-Pi-Seminara gauge 
theory [H]. The coupling of the BEC to its current 
gives rise to a number of exotic types of dynamics. We 
have shown how the presence of topological states cor- 
responding to persistent currents in a ring geometry de- 
pend on the number of particles. Also soliton solutions 
can be identified which are chiral in nature. It is cer- 
tainly tempting to draw analogies between the atomic 
system considered here and models of field theories de- 
scribing the fundamental forces between elementary par- 
ticles. From a quantum simulator point of view, perhaps 
the most intriguing aspect would be a generalisation of 
the mechanisms discussed here to a pseudo-spin situa- 
tion which can also support non-Abelian gauge poten- 
tials jSHl EI] • This will open up the possibility to address 
concepts from the standard model of particle physics and 
beyond, at an unprecedented level of sophistication. 
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